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Introduction 

One  critical  area  of  military  technology  that  is  relevant  to  both  the  Air  Force  and  the 
Strategic  Defense  Initiative  (SDI)  is  the  requirement  to  locate  ballistic  missiles,  particularly 
during  the  boost  phase  of  the  trajectory.  To  this  end,  codes  that  model  radiation  emitted 
from  rocket  plumes  have  been  developed  over  a  number  of  years  extending  back  to  the 
sixties.  With  the  advent  of  SDI,  the  requirements  on  these  codes  have  been  made  more 
stringent,  and  the  precision  required  in  their  predictions  has  increased,  in  large  part  due  to 
the  need  to  pinpoint  the  hard  body  within  the  extended  radiation  field  of  the  plume. 

The  need  for  kinetic  data  to  describe  high-altitude  phenomena  leading  to  plume 
signatures  has  been  evident  for  some  time.1  Since  the  primary  atmospheric  components  at 
altitudes  of  interest  are  N2  and  0(3P),  cross  sections  involving  these  species  with  exhaust 
molecules  are  of  particular  importance.  However,  experimental  results  for  processes 
involving  these  species  have  been  very  difficult  to  obtain.2  Because  experimental  data  for  a 
great  many  of  the  processes  of  interest  at  high  altitudes  are  simply  lacking,  reliable 
theoretical  data  have  become  essential  for  the  predictive  capabilities  of  plume  modelling 
codes  such  as  CHARM  and  SPURC. 

Theoretical  methods  are  applicable  to  a  wide  variety  of  molecular  systems  over 
large  velocity  ranges.  Thus  they  provide  a  means  to  obtain  much  of  the  required  data; 
primarily  excitation  cross  sections  and  rates  for  vibrational,  rotational  and  electronic 
energy  transfer  processes  (T-VRE)  and  for  chemical  reactions.  The  purpose  of  this 
research  was  to  develop  new  semiclassical  methods  to  improve  the  ability  of  theoretical 
methods  to  predict  kinetic  data  of  sufficient  accuracy  to  be  of  use  in  high-altitude  plume 
technology.  Virtually  all  of  the  theoretical  results  for  excitation  (T-VR)  processes  of 
interest  in  high  altitude  plume  kinetics  (i.e.  vibrational  and  rotational  excitation  of 
molecules  by  oxygen  atoms)  have  been  provided  by  the  Chemical  Dynamics  group  and 
their  collaborators  from  research  dating  back  to  the  late  seventies.3*4  This  research  was 
carried  out  under  sponsorship  of  the  Air  Force  Astronautics  Laboratory  (AFAL  - 
previously  AFRPL).  The  results  covered  a  number  of  systems,  including  H2O,  CO2,  HF, 
HC1,  and  CO.  The  theoretical  results  significantly  extended  the  data  base  of  available 
vibrational  excitation  cross  sections  and  provided  a  fairly  reliable  picture  of  the  relative 
magnitudes  of  the  cross  sections.  However,  of  necessity  at  that  time,  classical  mechanics 
was  used  for  the  dynamical  cross  section  calculations.  The  need  for  methods  to  more 
accurately  treat  vibrational  energy  transfer  processes  was  evident,  and  with  the  recent 
interest  in  UV  signature  phenomena,  a  method  applicable  to  electronic  excitation  was  also 
needed. 
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Successful  prediction  and  interpretation  of  UV-signatures  requires  detailed 
microscopic  information  about  the  responsible  collisional  processes.  One  of  the  central 
processes  in  this  connection  is  collision-induced  transitions  between  different  electronic 
states.  The  standard  quantum-mechanical  description  of  collisional  processes  is  based  on 
the  adiabatic  separation  of  nuclear  and  electronic  motions,  assuming  that  the  electrons 
instantaneously  adjust  to  nuclear  motions  because  of  the  large  difference  in  their  masses. 
Encounters  taking  place  on  a  single  adiabatic  surface  can  be  described  in  the  zeroth-order 
approximation  simply  by  the  classical  equations  of  motion.  However  transitions  from  one 
electronic  surface  to  another  are  essentially  quantum-mechanical  processes  which  have  no 
analogs  in  classical  mechanics.  The  task  of  numerically  solving  the  Schroedinger  equation 
governing  the  collisional  processes  in  question  is  still  an  unsolvable  computational  problem 
for  realistic  molecular  systems. 

Semiclassical  theory  fills  the  gap  between  the  quantum  and  classical  treatments. 
Approximating  the  quantum  mechanical  wavefunction  by  the  first  term  of  an  asymptotic 
expansion  leads  to  ordinary  first-order  differential  equations  for  the  semiclassical  solutions 
rather  than  the  more  difficult  partial  differential  equations  encountered  in  quantum 
mechanics.  This  approximation  is  valid  except  for  the  lowest  energies  and  is  sufficiently 
accurate  to  reproduce  specific  quantum-mechanical  features  which  are  completely  absent  in 
the  classical  picture.  In  particular  the  semiclassical  theory  is  capable  of  describing  the 
electronically  nonadiabatic  transitions  which  are  the  main  concern  in  the  present  research. 

There  is  however  an  essential  obstacle  that  has  been  hindering  progress  in  this  area 
for  many  years.  At  classical  turning  points  the  kinetic  energy  goes  to  zero  and  the 
semiclassical  solutions  are  singular.  In  one-dimensional  problems  this  situation  is 
corrected  using  connection  formulae.5  that  give  the  recipe  of  how  to  continue  the 
semiclassical  wavefunction  after  its  reflection  at  the  turning  point.  This  correction  exploits 
conservation  of  the  flux  to  connect  the  incoming  and  outgoing  solutions  but  is  insufficient 
in  the  multi-dimensional  case  as  there  may  be  exchange  of  flux  between  different  degrees  of 
freedom.  The  current  state-of-the-art  theory  in  this  area  is  mainly  developed  at  the  level  of 
formal  theorems  which  prove  existence  of  solutions,  but  do  not  provide  any  practical 
numerical  algorithms  to  find  them. 


Objectives 

The  objectives  of  Phase  I  were  (1)  to  identify  relevant  electronic  structure  input  for 
UV -signature-producing  collisional  processes  involving  excited  NO,  OH  and  N2 
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molecules,  and  (2)  to  develop  semiclassical  techniques  capable  to  provide  the  necessary 
state-to-state  dynamical  information  for  collisional  processes  relevant  to  SPURC. 

Task  1  involves  examining  the  molecular  correlation  diagrams  for  collisional 
processes  between  relevant  atmospheric/plume  species  in  order  to  identify  the  associated 
quantum  chemical  potential  energy  surfaces  and  electronic  couplings,  operative  in  the 
energy  range  of  interest  Task  2  involves  theoretical  developments  that  enable  semiclassical 
wavefunction  computations  for  multidimensional  electronically  inelastic  systems. 

Status  of  Research 

Significant  progress  on  both  the  tasks  of  Phase  I  is  reported  by  presenting  detailed 
results.  Most  notably,  task  2  has  resulted  in  a  significant  progress  in  overcoming  two  most 
serious  problems  of  the  nonadiabatic  semiclassical  theory,  making  possible  its  practical 
implementation  for  the  quantitative  description  of  processes  leading  to  UV  radiation.  In 
addition,  an  improved  method  of  treating  T-V  excitation  on  a  single  adiabatic  surface  was 
also  developed. 

The  adiabatic  separation  of  nuclear  motions  in  the  semiclassical  picture  made  it 
possible  to  connect  incoming  and  outgoing  semiclassical  solutions,  reproducing  correctly 
interference  effects  between  different  trajectories.  Validity  of  the  adiabatic  separation  of 
streamline  and  vibrational  motions  was  demonstrated  by  comparing  the  semiclassical  and 
exact  quantum-mechanical  results  for  vibrational  energy  transfer  in  an  atom-diatom 
collision  in  the  IOS  approximation6. 

Including  the  adiabatic  hypothesis  into  the  Self-Consistent  Eikonal  Method 
(SCEM)7  gave  us  a  practical  tool  for  treating  electronically  inelastic  processes  with  a  few 
degree  of  freedom.  As  another  test  of  applicability  of  the  single  average  potential  (the  so- 
called  Ehrenfest  potential  obtained  by  weighting  the  adiabatic  potential  functions  with  the 
probabilities  for  the  system  to  be  in  the  appropriate  electronic  state)  we  studied  vibrational 
excitation  of  the  diatomic  molecule  in  an  atom-diatom  collision  in  the  IOS  approximation  by 
expanding  the  Schroedinger  equation  in  the  vibrational  basis  set. 

Phase  I  also  resulted  in  the  development  of  a  new  extended  approach  called 
SCACEM  (Self-Consistent  Adiabatically-Corrected  Eikonal  Method)  for  a  more  accurate 
treatment  of  electronic  inelasticity  in  molecular  collisions;  this  goes  beyond  the  Ehrenfest 
effective  potential  approximation  employed  in  SCEM  without  compromising  on  the 
advantages  of  the  latter  method.  The  following  executive  summary  emphasizes  the  relevant 
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implications  of  our  computational  approach  for  future  SPURC-related  investigations. 
Phase  I  results  are  overviewed  in  a  following  section.  Two  manuscripts  (one  on  the 
extended  formal  theory  and  one  on  numerical  studies)  are  included  as  appendices  in  this 
report. 

Executive  Summary 

An  analysis  of  the  correlation  diagrams  for  dominant  collisional  pathways  leading  to 
excited  states  of  the  three  candidate  molecules  (NO,  OH,  and  N2)  was  conducted  to  identify 
a  feasible  system  for  quantum  chemistry  input.  The  problem  of  OH(2L+)  radical  formation 
in  the  reactive  process: 

0(3P)  +  H2  -»  OH(2E+)  +  H(2S)  (A) 

emerges  as  an  immediately  feasible  case  for  quantum  chemical  characterization  of  potential 
energy  surfaces  and  interelectronic  state  couplings.  Clearly,  the  dynamical  treatment  of 
process  (A)  requires  a  computational  description  of  electronic  nonadiabaticity  as  well  as 
reaction  (nuclear  rearrangement)  leading  to  an  electronically  excited  OH(2E+)  product 
during  collisions. 

In  Phase  I  we  successfully  overcame  the  most  serious  principal  defect  of  the 
semiclassical  theory  caused  by  its  break-down  at  classical  turning  points.  We  developed  a 
new  technique  to  construct  semiclassical  wavefunctions  for  realistic 
potential  surfaces  regardless  of  the  complexity  of  the  latter.  For  electronically 
inelastic  processes  this  new  technique  is  included  as  an  integral  part  either  in  the  SCEM 
(whose  effectiveness  has  been  proved  earlier  and  is  additionally  supported  by  new 
calculations  performed  during  Phase  I)  or  in  the  aforementioned  new  theory,  SCACEM, 
which  is  expected  to  be  applicable  even  in  those  cases  when  potential  surfaces  in  different 
electronic  states  significantly  differ  from  each  other. 

Mathematical  arguments  presented  in  support  of  the  new  formalism  make  it 
especially  attractive  for  the  semiclassical  description  of  reactive  encounters.  In  fact, 
although  the  latter  may  be  formally  treated  within  SCEM,  they  may  involve  large  relative 
deviations  of  different  potential  energy  surfaces,  thereby  leading  to  errors  from  employing 
a  simple  average  Ehrenfest  potential.  The  improved  theory,  SCACEM,  generated  in  Phase 
I  as  a  feasible  strategy  to  accurately  treat  processes  involving  both  electronic  inelasticity  and 
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reaction,  overcomes  this  difficulty  by  running  trajectories  separately  on  each  potential 
surface. 

Another  complication  associated  with  treatment  of  reactive  encounters  by  means  of 
a  single  effective  potential  comes  from  the  fact  that  stretching  motions  of  interatomic  bonds 
undergo  a  sharp  change  when  bonds  break;  in  order  to  select  coordinates  of  quasi-periodic 
motions  one  has  to  refer  the  vibrations  to  the  appropriate  reaction  path.  Since  the  reaction 
path  varies  from  one  potential  surface  to  another,  it  requires  use  of  different  variables  for 
different  potential  surfaces;  this  is  relatively  easy  to  implement  in  the  frame  of  the  new 
theory.  The  new  theory  also  enables  to  describe  tunneling  effects  along  the  streamline 
coordinate,  including  a  certain  type  of  resonances. 

During  Phase  I  the  applicability  of  the  developed  semiclassical  theory  has  been 
tested  and  validated  for  a  two-dimensional  representation  of  vibration?!  excitation  in  O+HF 
and  Na+H2  collisions  for  which  we  can  generate  exact  quantum  results.  These  successful 

applications  encourage  further  exploration  of  the  technique  in  realistic  multidimensional 
applications. 

In  conclusion,  Phase  I  research  has  (1)  revealed  that  reliable  ab  initio  quantum 
chemical  input  can  be  obtained  for  SPURC-related  systems,  (2)  established  that  elementary 
computational  modules  required  to  build  the  necessary  dynamical  treatments  for  SPURC- 
related  rate  computations  are  reliable  (Appendix  2),  (3)  resulted  in  a  significant  theoretical 
break-trough  in  development  of  the  methodology  capable  to  provide  the  nuclear  state-to- 
state  information  that  both  SPURC  and  CHARM  ultimately  need.  It  is  also  clear  that 
further  systematic  method  refinement  is  essential  since  this  type  of  computational  input  is 
valuable  to  the  USAF  due  to  the  cost  and  difficulties  for  experimental  determination  of  the 
required  data.  Towards  this  end  an  extended  semiclassical  theoretical  methodology  based 
on  the  eikonal  wavefunction  ansatz  (Appendix  1)  should  be  employed. 

Phase  I  Results 

1 .  SPURC-related  system  survey 

The  three  candidate  molecules  considered  were  OH,  NO,  and  N2.  Correlation 
diagrams  for  the  following  processes  leading  to  the  production  of  excited  species  were 
examined: 

0(3P)+  N2  NO(2E+)  +  N(4S) 
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0(3P)  +  H2  -4  OH(2L+)  +  H(2S) 

0(3P)+  N2(X1Zg+)  — »  N2(1Lu+)  +  0(3P) 


Details  of  the  pertinent  electronic  states,  dominant  interelectronic  state  couplings,  and  the 
likely  dynamical  pathways  for  the  production  of  the  excited  species  were  analyzed.  The 
feasibility  of  obtaining  the  quantum  chemical  input  was  examined,  and  the  conclusion  is 
that  accurate  surfaces  can  be  obtained  for  the  (O.H2)  system,  with  useable  although  less 

accurate  surfaces  for  other  systems.  We  also  point  out  that  some  of  the  processes  emerging 
from  the  recent  kinetics  work8  at  Aerospace  Corporation  reveal  important  electronic  energy 
transfer  processes  relevant  to  UV  signatures,  and  we  believe  that  a  theoretical  treatment  can 
provide  useful  information  to  augment  the  experimental  data.  Figure  1  is  a  schematic  of  the 
correlation  of  fragment  molecule  states  important  for  the  (O,^)  system.  Figure  2  provides 
similar  information  for  the  (O,!^)  system. 

2.  Eikonal  Wavefunction  Method 
Theory  Development 

The  formal  issues  to  be  examined  in  Phase  I  were: 

(i)  development  of  prescriptions  for  continuation  of  global  semiclassical 
wavefunctions  through  classical  turning  points, 

(ii)  development  of  a  practical  numerical  algorithm  to  extract  final  populations  of 
vibrational  levels  from  the  semiclassical  outgoing  solution, 

(iii)  implementation  of  the  single-surface  developments  (i),  (ii)  in  the  multichannel 
theory  with  valid  use  of  the  Ehrenfest  effective  potential  (SCEM), 

(iiii)  exploration  of  means  to  go  beyond  the  assumptions  underlying  SCEM  to  take 
into  account  specific  topology  of  the  adiabatic  potential  surfaces,  especially  responsible  for 
tunneling  and  related  resonances  which  gave  rise  to  the  new  multichannel  semiclassical 
theory  referred  to  as  SCACEM. 


The  details  of  Phase  I  formal  developments  on  all  these  issues  are  covered  in  the  formalism 
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paper  included  as  Appendix  1.  Here  we  briefly  summarize  the  highlights  of  these 
developments  in  less  mathematical  terms. 


Figure  1.  Electronic  correlation  diagram  (Cs  symmetry)  for  low-lying  states  of  the 
OC^P)  +  H2  system. 


Adiabatic  Approximation  of  Nuclear  Motions  in  the  Semiclassical  Picture 

When  the  primary  reason  for  failure  of  classical  mechanics  is  due  to  interference 
effects  between  amplitudes  along  different  trajectories  (this  is  typically  the  case  when 
tunneling  and  curve  crossing  are  unimportant),  semiclassical  techniques  have  shown  the 
potential  to  provide  a  very  efficient  computational  route  that  is  more  accurate  than  classical 
mechanics  and  less  expensive  than  quantal  methods.  Their  main  advantage  is  that 
propagation  of  wavefunctions  along  classical  trajectories  is  governed  by  a  set  of  ordinary 
differential  equations  of  the  first  order  which  is  much  easier  to  solve,  compared  with  the 
Schroedinger  equation.  There  is  however  a  serious  obstacle  that  damped  any  progress  in 
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Figure  2.  Electronic  correlation  diagram  (Cs  symmetry)  for  states  of  the  (O(^P), 
N2}  system. 


this  direction  for  many  years  -  the  semiclassical  approximation  breaks  down  for 
low  velocities  such  as  those  experienced  when  the  system  hits  a  potential  wall  or  barrier. 
In  one-dimensional  systems  the  problem  is  easily  corrected  using  standard  WKB 
connection  formulae  to  carry  the  solutions  through  those  turning  points  into  classically 
forbidden  regions.  Multidimensional  systems  present  a  much  more  difficult  problem 
because  the  probability  density  may  flow  from  one  degree  of  freedom  to  another  and  hence 
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one  can  no  longer  take  advantage  of  flux  conservation  to  connect  incoming  and  outgoing 
solutions. 

At  the  present  there  are  no  practical  numerical  algorithms  for  the  general  solution 
of  this  problem.  To  by-pass  the  difficulties  coming  from  the  turning  points  we  restrict  our 
analysis  to  adiabadcally  separable  motions  so  that  one  can  apply  the  WKB-connection 
formulae  for  each  degree  of  freedom,  namely,  we  subtract  tc/2  from  the  Hamilton 
characteristic  function  each  time  when  one  of  the  conjugate  momenta  changes  its  sign. 

As  discussed  in  Appendix  1  there  are  two  practical  problems  which  can  be 
immediately  treated  in  such  a  way:  nonreactive  atom-diatom  scattering  and  scattering  of 
diatomic  molecules  by  a  surface.  To  illustrate  the  specifics  of  adiabatic  separation  of 
nuclear  motions  in  the  semiclassical  theory,  let  us  consider  atom-diatom  collisions  with 
zero  angular  momentum.  Relative  oscillations  of  two  atoms  and  the  nutation  of  the  diatom 
give  us  two  quasi-periodic  degrees  of  freedom.  As  stretching  vibration,  described  by  the 
variable  r,  usually  has  a  much  larger  frequency  than  the  nutation,  we  treat  the  former  as  a 
quasi-periodic  motion  with  a  frequency  parametrically  dependent  on  the  nutadonal  angle  0 
and  the  streamline  coordinate  R.  On  the  other  hand,  the  nutadon  is  expected  to  be  faster 
than  streamline  motion  so  that  we  can  describe  it  as  quasi-periodic  motion  with  the 
frequency  dependent  parametrically  on  R.  We  thus  neglect  the  effect  of  small-amplitude 
oscillations  of  the  interatomic  distance  r  on  the  nutation  of  the  diatomic  and  treat  the 
diatomic  as  a  rigid  linear  rotator  in  this  particular  context.  It  is  important  that  one  needs 
these  assumptions  only  to  formulate  the  connection  formulas  for  the  asymptotic  solutions 
and  we  do  not  make  any  approximations  in  the  classical  hamiltonian  itself  to  be  consistent 
with  these  assumptions. 

An  extension  of  the  developed  technique  to  describe  reactive  scattering  of  two 
polyatomic  fragments  is  a  more  sophisticated  problem  which  will  be  treated  within  the 
reaction-path  formalism. 

Numerical  Evaluation  of  the  S-Matrix  Elements 

Singularity  of  semiclassical  wavefunctions  in  turning  points  gives  rise  to  another 
problem  —  a  very  slow  convergence  of  the  their  overlaps  with  quantum  mechanical 
vibrational  wavefunctions  used  to  evaluate  final  populations  of  vibrational  levels.  To 
overcome  the  numerical  difficulties  associated  with  it  we  developed  a  new  mathematical 
scheme  by  formally  transforming  the  (N-l)-dimensional  integral,  used  for  evaluation  of  the 
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S-matrix  element  for  a  scattering  process  with  N-l  quasi-periodic  degrees  of  freedom,  into 
a  N-dimensional  integral  such  that  singularities  of  the  semiclassical  wavefunction  is 
compensated  by  zeros  of  the  Jacobian  in  the  N-dimensional  volume  element.  The 
appropriate  mathematical  details  are  given  in  Section  VI  in  Appendix  1. 

Self-Consistent  Eikonal  Method  (SCEM)  in  Applications  to  Multidimensional 

Problems 

Finally,  for  electronically  inelastic  systems,  we  recommend  that  the  specified 
multidimensional  nuclear  wavefunction  be  combined  with  the  self-consistently  calculated 
electronic  amplitudes  to  obtain  the  overall  (electronic+nuclear)  state-to-state  transition 
amplitudes  including  interference  effects.  Section  IV  in  Appendix  1  outlines  the  main 
features  of  the  SCEM  propagating  trajectories  in  the  field  of  a  single  average  potential.  It 
was  shown  that  SCEM  corresponds  to  a  very  naturai  physical  picture,  namely,  the  flows  of 
the  probability  density  on  different  adiabatic  potential  surfaces  are  approximated  by  the 
average  flow  obtained  by  weighting  the  individual  flows  with  the  probabilities  for  the 
system  to  be  in  the  appropriate  electronic  state.  We  also  revealed  some  additional  terms  in 
the  potential  forces  governing  nuclear  motions  in  the  case  of  a  system  with  two  or  more 
degrees  of  freedom  and  suggested  the  numerical  algorithm  for  their  evaluation  (see 
Appendix  A  in  the  enclosed  formalism  paper  referred  to  above  as  "Appendix  1"). 
Significance  of  those  terms  for  SPURC-related  applications  is  expected  to  be  studied  in 
more  details  5n  future  research. 

SCEM  may  be  employed  for  the  nonreactive  curve-crossing  parts  to  propagate  the 
wavefunction  amplitudes  up  to  a  reactive  region  where  the  single-surface  methods  such  as 
variational  transition  state  theory  (VTST)  can  continue  the  computations  to  yield  the 
required  rate  constants.  An  advantage  of  using  VTST  is  that  one  avoids  the  necessary 
developments  to  treat  the  tunneling  in  directions  perpendicular  to  the  reaction  path. 
Remember  that  VTST  has  developed  the  reliable  recipes  to  treat  this  effect  with  a  sufficient 
accuracy  whereas  including  it  in  the  semiclassical  technique  is  a  relatively  complicated 
proMem.  Ideally,  for  SPURC  problems,  it  will  be  necessary  to  propagate  SCEM 
amplitudes  through  the  nonadiabatic  region  and  to  use  the  constructed  semiclassical 
wavefunctions  to  extract  initial  populations  of  states  for  the  single-channel  problem  in  the 
reactive  region. 
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Self  Consistent  Adiabatically-Corrected  Eikonal  Method  (SCACEM) 

The  necessary  state-to  state  information  can  be  obtained  from  the  new  extended 
theory,  SCACEM,  which  exploits  the  eikonal  ansatz  for  nuclear  wavefunctions  in  each 
electronic  state,  where  each  wavefunction  is  propagated  now  along  its  own  trajectory  run 
on  the  appropriate  adiabatic  potential  surface.  The  new  formalism  is  the  most  significant 
development  arising  from  Phase  I  for  SPURC  and  was  developed  to  address  situations 
where  the  Ehrenfest  effective  potential  used  in  SCEM  can  lead  to  inaccuracies.  For 
example,  it  may  be  anticipated  that  the  Ehrenfest  effective  potential  employed  in  SCEM  can 
be  inaccurate  when  a  reactive  barrier  region  appears  only  on  one  of  the  electronic  surfaces. 

One  of  the  advantages  of  propagating  classical  trajectories  on  the  particular  adiabatic 
surfaces  as  prescribed  by  SCACEM  is  that  additional  approximations  are  made  only  in  the 
nonadiabatic  region.  The  equations  beyond  this  region  are  the  exact  limit  of  the 
Schroedinger  equation  in  the  adiabatic  representation  as  h  tends  to  zero  and  hence  their 
solutions  reproduce  asymptotically  specific  features  of  the  exact  wavefunctions  caused  by 
local  topology  of  the  adiabatic  potential  surfaces.  In  particular,  the  new  technique  makes  it 
possible  to  trei  t  tunneling  along  the  streamline  coordinate  R  in  electronically  inelastic 
processes  provided  it  takes  place  beyond  the  region  of  strong  nonadiabatic  couplings. 

Let  us  consider  as  an  example  three  lowest  adiabatic  surfaces  for  Na+H^  C2V 

encounters.  Fig.3  represents  the  potential  surfaces  reported  by  Eakes9  for  the  H-H  distance 
equal  to  1.4aQ.  We  consider  collisions  with  energies  2.2-2.5eV  between  the  bottom  and  the 

barrier  on  the  upper  curve.  We  start  propagation  at  a  large  value  of  R  where  two  excited 

electronic  states  strongly  interact  with  each  other  so  that  the  probability  density  flows  from 

the  initially  populated  excited  state  to  the  second  one.  In  the  region  where  the  trajectory  run 

on  the  upper  surface  hits  the  potential  barrier  the  nonadiabatic  interactions  become 

negligible  and  we  can  evaluate  the  reflection  coefficient  by  using  the  standard  formulae  of 

the  one-dimensional  tunneling  theory10.  As  for  the  trajectory  on  the  lower  excited  electronic 
surface  it  comes  back  after  hitting  the  potential  wall  at  R=3aQ  and  then  both  flows  start 

interfering  again  with  each  other.  Obviously  if  the  upper  potential  well  has  quasi-bound 
states  they  manifest  as  resonances  in  the  scattering  amplitudes  evaluated  according  to  the 
outlined  recipe. 
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Figure  3.  Na  +  H2  potentials  for  C2V  approach  of  Na. 

Computational  implementation  of  SCACEM  and  tests  to  reveal  its  potential  for  the 
problems  on  hand  should  be  considered  as  one  of  the  most  important  tasks  in  the  future 
research. 

Numerical  Work 

The  best  route  to  establish  the  feasibility  of  approximate  semiclassical  wavefunction 
computations  is  by  direct  verification  in  applications  to  specific  benchmark  systems  for 
which  exact  quantum-mechanical  results  are  available.  To  test  the  developed  algorithms  we 
used  the  two-dimensional  model  obtained  from  the  multidimensional  Schroedinger  equation 
in  the  IOS  approximation.  The  exact  calculations  with  such  a  model  for  vibrational 
excitation  of  HF  in  collisions  with  oxygen  have  been  reported  earlier  by  us  and  moreover 
any  new  necessary  results  could  be  easily  generated  by  using  the  available  quantum  code. 
For  a  similar  reason  we  chose  to  study  the  Na+H2  as  another  example.  This  system  is  more 

relevant  to  the  SPURC-related  programs  as  the  hydrogen  molecules  may  undergo 
nonadiabatic  transitions.  However  the  latter  feature  makes  exact  quantum  mechanical 
calculations  much  more  complicated  even  in  their  simplest  two-dimensional  reduction  so 
that  the  converged  exact  solutions  could  be  obtained  up  to  now  only  for  vibrational 
excitation  of  H2  in  the  ground  electronic  state.  We  thus  had  to  restrict  our  analysis  at  this 

stage  only  to  those  electronically-elastic  processes.  It  is  worth  emphasizing  that,  despite 
the  simplicity  of  the  model  used  for  our  tests  it  has  inherited  the  main  problem  of  the 
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multidimensional  semiclassical  theory  —  its  break-down  on  the  boundary  of  the  classically 
allowed  region. 

Compared  with  the  one-dimensional  case,  the  probability  density  is  now 
redistributed  between  different  trajectories,  leading  to  nonconservation  of  the  flux  at  the 
moment  the  trajectory  hits  the  boundary  surface.  As  an  extension  of  our  approach  to  more 
realistic  systems  requires  solution  of  mathematical  difficulties  similar  to  those  we  have 
addressed  in  this  research,  the  method  does  have  the  potential  to  provide  useful  information 
about  the  dynamics  of  real  collisional  processes.  The  accuracy  of  this  information  for 
realistic  multidimensional  applications  is  now  being  tested  for  rotationally  inelastic 
scattering  of  NO  by  an  uncorrugated  silver  surface  —  the  example  earlier  treated  by  us  using 
rotonic  expansions. 

At  the  very  outset,  the  key  feasibility  issue  to  be  addressed  in  Phase  I  was 
identified  to  be  the  validation  of  the  semiclassical  eikonal  wavefunction  route  for  computing 
nuclear  state-to-state  collisional  excitation  propensities.  Two  examples  of  collisional 
excitation  of  vibrational  excited  states  were  chosen  for  which  exact  quantum  limit  results 
were  available.  These  are  O+HF  and  Na+H2  collisions  treated  within  a  rotational  infinite 

order  sudden  (IOS)  model.  The  QIOS  (quantal  infinite  order  sudden)  models  are  typically 
the  largest  available  computations  for  most  systems  owing  to  the  proliferation  of  rotational 
channels  and  they  do  contain  the  essential  physics  of  vibrational  excitation  dynamics  which 
is  what  the  eikonal  wavefunction  was  being  tested  for  its  ability  to  describe. 

Phase  I  numerical  work  focused  on  two  distinct  relevant  subtasks: 

(1)  method  development  for  computing  S-matrices  using  semiclassical  eikonal 
wavefunction  solutions  constructed  from  nuclear  trajectories,  thus  including  quantal 
interference  effects  for  nuclear  motions  (numerical  tests  were  based  on  an  implementation 
for  O+HF  collisions  leading  to  vibrational  excitations),  and 

(2)  numerical  tests  on  a  SCEM  treatment  of  atom-diatom  collisions  which  used  a 
vibronic  expansion  (thus  quantum  mechanically  treating  vibrations)  while  classically 
treating  translation  and  within  the  rotational  infinite  order  sudden  model  of  0+HF(v) 
0+HF(v')  as  well  as  Na+H2(v)  -»  Na+H2(v')  collisions.  The  treatment  employed  here  is 

often  known  as  the  coupled  states  treatment  and  is  another  way  of  obtaining  nuclear  state- 
to-state  information  along  with  electronic  state-to-state  information. 
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The  dynamical  model  (1)  is  called  the  SIOS  (for  semiclassical  infinite-order- 
sudden)  model  and  (2)  is  called  the  RSVE  (for  rotational  sudden  with  vibronic  expansion) 
model.  Appendix  2  has  a  preliminary  draft  for  a  paper  resulting  from  these  numerical 
studies.  The  key  result,  represented  by  the  data  in  Fig.  4,  is  that  even  a  modest  calculation 
(employing  20  trajectories)  within  the  SIOS  model  is  able  to  obtain  very  good  prediction  of 
microscopic  vibrational  excitation  probabilities,  thereby  implying  that  reliable  state-to-state 
vibrational  excitation  cross  sections  will  result  from  such  a  computation. 


O+HF  Vibrational  Excitation 


Figure  4.  Benchmark  based  on  rotational  infinite  order  sudden  (IOS)  angle-dependent 

probabilities  for  v=0  ->  v=l  vibrational  excitation  in  O(^P)+HF(v=0,j=0)  collisions  at 
Etr=3.0eV,  for  total  angular  momentum,  L=0.  QIOS  denotes  exact  quantum  IOS 

calculation,  RSVE  denotes  rotational  sudden  with  vibrational  expansion,  and  SIOS  denotes 
semiclassical  infinite  order  sudden.  SIOS  results  are  based  on  semiclassical  eikonal 
wavefunctions  (obtained  here  from  a  meagre  20  classical  trajectories  and  hence  expected  to 
improve)  and  promise  the  practically  most  attractive  route  for  these  excitation  propensities 
in  signature  problems. 


The  results  shown  in  Fig.  4  are  for  relatively  low  values  of  probabilities  (near 
threshold)  and  this  is  a  crucial  aspect  required  of  a  dynamical  method  to  be  useful  for 
signature  problems:  the  rationale  is  that  larger  Boltzmann  weights  come  in  for  lower 
energies  and  thus  the  low-intermediate  energy  contribution  has  a  significant  role  to  play  in 
the  observed  IR/UV  signature.  Such  velocity  components  are  present  to  some  degree  in 
plume  flowfields  involving  even  hypersonic  vehicles.  The  present  method  offers  a 
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significant  improvement  over  the  previously  available  classical  trajectory  method  in  this 
regard.  The  RSVE  method  is  also  seen  to  be  successful,  but  here  we  stress  the  new  SIOS 
route  which  is  much  more  economical  and  has  the  potential  to  scale  up  to  much  larger 
polyatomic  problems,  a  regime  where  the  basis  set  expansion  methods  become  unwieldy. 
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Appendix  1 


SEMICLASSICAL  METHODOLOGY  FOR  COMPUTING  MULTICHANNEL  EIKONAL 
WAVEFUNCTIONS  IN  MOLECULAR  COLLISIONS:  A  REFORMULATION  AND 

EXTENSION 


by 

G.  A.  Natanson  and  P.  K.  Swaminathan 
Chemical  Dynamics  Corporation 
Upper  Marlboro,  MD  20772 


I.  Introduction 


The  standard  quantum-mechanical  description  of  collisional  processes  is  based  on  the 
adiabatic  separadon  of  nuclear  and  electronic  modons,  assuming  that  the  electrons  instantaneously 
adjust  to  neaclear  modons  because  of  a  large  difference  in  their  masses.  Encounters  taking  place  on 
a  single  adiabadc  surface  can  be  described  in  the  zeroth  approximation  simply  by  the  classical 
equations  of  motion.  However  transitions  from  one  electronic  surface  to  another  are  essentially 
quantum-mechanical  processes  which  have  no  analogs  in  classical  mechanics.  As  numerical 
solution  of  the  Schroedinger  equation  is  near  the  capability  limit  of  modem  computers  even  for 
four-atom  collisions  on  a  single  potential  surface1  the  prospects  of  obtaining  the  reliable 
information  for  nonadiabatic  processes  do  not  seem  rather  promising  at  the  present. 

The  gap  between  the  quantum  and  classical  treatments  can  be  however  filled  by  the 
semiclassical  theory.  Construction  of  semiclassical  solutions  is  simpler  compared  with  solving  the 
Schroedinger  equation  because  they  are  governed  by  ordinary  differential  equations  of  the  first 
order  (instead  of  partial  differential  equations  exploited  in  quantum  mechanics).  On  the  other  hand, 
the  semiclassical  sheory  is  accurate  enough  to  reproduce  specific  quantum-mechanical  features 
which  are  completely  absent  in  the  classical  picture.  In  particular  the  semiclassical  theory  turns  out 
to  be  powerful  enough  to  describe  the  nonadiabatic  transitions  which  are  the  main  concern  in  the 
present  research. 

There  is  however  an  essential  obstacle  that  has  been  damping  a  progress  in  this  area  for 
many  years.  The  semiclassical  equations  imply  that  higher-order  corrections  are  negligible 
compared  with  the  kinetic  energy  the  system  has  so  that  they  break  down  each  time  when  the  latter 
becomes  small  as  it  takes  place  at  classical  turning  points  where  semiclassical  wavefunctions  are 
singular.  The  one-dimensional  problem  has  the  well-known  solution  in  terms  of  the  so-called 
WKB  method2  that  gives  the  recipe  of  how  to  continue  the  semiclassical  wavefunction  after  its 
reflection  in  the  turning  point.  This  recipe  however  essentially  exploits  conservation  of  the  flux  to 
connect  the  incoming  and  outgoing  waves  which  is  insufficient  in  the  multi-dimensional  case  as 
there  may  be  exchange  of  the  probability  density  between  different  degrees  of  freedom.  The 
current  state-of-the-art  theory  in  this  area  is  mainly  developed  at  the  level  of  formal  theorems3 
which  prove  existence  of  solutions,  but  do  not  provide  any  practical  numerical  algorithms  to  find 
them.  We  can  cite  only  a  few  works4’5  where  the  semiclassical  scattering  theory  in  many 
dimensions  was  brought  to  the  particular  applications. 

An  additional  problem  arises  when  one  applies  the  semiclassical  technique  to  nonadiabatic 
processes.  It  comes  from  the  interference  between  wavefunctions  propagated  on  different  potential 
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surfaces.  To  take  into  account  this  effect  one  needs  to  locate  on  each  surface  the  trajectory  that  goes 
through  the  given  point  in  space  which  is  numerically  unaccessable. 

Up  to  now  a  much  more  attention  has  been  to  given  an  alternative  direction  exploiting 
angle-action  variables,  in  following  pioneering  works  of  Marcus6  and  Miller7  however  the  reported 
success  in  this  direction  is  also  very  limited.  Even  disregading  the  fact  that  this  most  popular 
approach  is  related  with  the  Schroedinger  equation  only  through  the  Correspondence  Principle  and 
hence  it  remains  unclear  to  what  extent  its  predictions  are  equivalent  to  those  of  the  asymptotic 
semiclassical  theory3’8  there  is  a  pure  practical  reason  for  development  of  the  semiclassical  theory 
in  space  —  representation  of  realistic  potential  surfaces  in  angle-action  variables  is  a  complicated 
computational  problem  that  makes  the  whole  scheme  unfeasible  for  applications. 

The  present  research  makes  a  significant  progress  in  overcoming  two  most  serious 
problems  of  the  nonadiabatic  semiclassical  theory,  making  possible  its  practical  implementation. 
The  assumption  about  adiabatic  separation  of  nuclear  motions  in  the  semiclassical  picture  allows  us 
to  connect  incoming  and  outgoing  semiclassical  solutions,  reproducing  correctly  interference 
effects  between  different  trajectories.  Validity  of  the  adiabatic  separation  of  streamline  and 
vibrational  motions  is  demonstrated  by  comparing  the  semiclassical  and  exact  quantum-mechanical 
results  for  the  two-dimensional  model  describing  vibrational  energy  transfer  in  an  atom-diatom 
collision  in  the  IOS  approximation9. 

Including  the  adiabatic  hypothesis  into  the  Self-Consistent  Eikonal  Method  (SCEM)10  gave 
us  a  practical  tool  for  treating  electronically  inelastic  processes  with  a  few  degree  of  freedom.  It 
was  shown  that  SCEM  corresponds  to  a  very  natural  physical  picture,  namely,  the  flows  of  the 
probability  density  on  different  adiabatic  potential  surfaces  are  approximated  by  the  average  flow 
obtained  by  weighting  the  individulal  flows  with  the  probabilities  for  the  system  to  be  in  the 
appropriate  electronic  state.  We  also  revealed  some  additional  terms  in  the  potential  forces 
governing  nuclear  motions  in  the  case  of  a  system  with  two  or  more  degrees  of  freedom  and 
suggested  the  numerical  algorithm  for  their  evaluation. 

We  also  developed  a  new  extended  approach  called  SCACEM  (Self-Consistent 
Adiabatically-Corrected  Eikonal  Method)  for  a  more  accurate  treatment  of  electronic  inelasticity  in 
molecular  collisions;  this  goes  beyond  the  Ehrenfest  effective  potential  approximation  employed  in 
SCEM  without  compromising  on  the  advantages  of  the  latter  method.  The  new  extended  theory 
also  exploits  the  eikonal  ansatz  for  nuclear  wavefunctions  in  each  electronic  state  but  each 
wavefunction  is  propagated  now  along  its  own  trajectory  run  on  the  appropriate  adiabatic  potential 
surface.  One  of  the  advantages  of  propagating  classical  trajectories  on  the  particular  adiabatic 
surfaces  as  prescribed  by  SCACEM  is  that  additional  approximations  are  made  only  in  the 
nonadiabatic  region.  The  equations  beyond  this  region  are  the  exact  limit  of  the  Schroedinger 
equation  in  the  adiabatic  representation  as  h  tends  to  zero  and  hence  their  solutions  reproduce 
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asymptotically  specific  features  of  the  exact  wavefunctions  caused  by  local  topology  of  the 
adiabatic  potential  surfaces.  In  particular,  the  developed  technique  makes  it  possible  to  treat 
tunneling  along  the  streamline  coordinate  R  in  ellectronically  inelastic  processes  provided  it  takes 
place  beyond  the  region  of  strong  nonadiabatic  couplings. 

II . Semiclassical  multichannel  approximation  in  the  mixed 
adlabatlc/diabatic  representation 


Let  us  consider  the  multichannel  problem  described  by  a  set  of  linear  partial  differential 
equations  of  the  second  order: 


(  -  ~Aq  -  e]i  +  V(Q)  -  h2  l (Q)  •  V  + 


H(Q)]  y  =  0  , 


(2.1) 


where  V(Q)  is  the  column  formed  by  the  functions  \|/j(Q),  \y2(Q) . Yn(Q)  sought  f°r. 


4,=  If1/2 


MQ  =  *  g 

VW 


dQ 


dQ 


(2.2) 


and  V  are  respectively  the  covariant  Laplacian  (see  next  Section  for  comments)  and  the  gradient 
in  the  space  of  curvilinear  coordinates  Q  (see  next  Section  for  comments),  V(Q)  is  a  n  x  n 

diagonal  matrix  having  some  potentials  Vy(Q)  as  its  diagonal  elements,  H(Q)  is  the  matrix  of 

potential  coupling  having  zeros  at  its  main  diagonal;  the  components  of  the  vector  F  are  n  x  n  real 
antisymmetric  matrices  describing  nonadiabatic  couplings.  In  the  adiabatic  representation  Vy(Q) 
are  the  adiabatic  potentials  and  the  matrix  H(Q)  s  0  so  that 


1  +  V(Q)  -  h 


2 


F(Q).F 


\j/  =  0  . 


(2.3) 


The  main  advantage  of  the  adiabatic  representation  is  that  the  last  term  in  brackets  in  (2.3)  vanishes 
in  the  limit  as  h  -*  0.  As  a  result  we  come  to  the  Hamilton-Jacobi  equations  governed  by  the 
adiabatic  potentials.  We  can  thus  treat  the  coupling  term  as  a  perturbation  .  We  show  below  that 
one  needs  an  assumption  of  such  a  kind  when  continuing  the  semiclassical  wavefunction  through 
turning  points  of  quasiperiodic  motions.  Another  alternative  is  a  system  with  only  potential 
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couplings  which  are  relatively  small  compared  with  the  potentials  Vy(Q).  In  the  one-dimensional 

case  it  is  sufficient  to  require  that  the  couplings  are  negligible  near  the  turning  point  of  the 
streamline  motion.  We  come  to  the  problem  of  this  type,  when  treating  quasiperiodic  degrees  of 
freedom  quantum-  mechanically  by  representing  the  Schroedinger  equation: 


f 


±i.  + 

2  ^Q  + 


V(Q)  -  E  V  =  0 


(2.4) 


in  the  matrix  form: 


ff  h2  d2  1 

' 

-  ,-E 

1  +  V(R)  +  H(R) 

.V  2  dR2  ) 

|=  _  - 

V  =  0  , 


(2.5) 


where  R  is  the  streamline  coordinate  which  is  assumed  to  be  orthogonal  to  all  others.  As  a  result 
Eq.(2.4)  takes  the  form: 


(  d2  1 

| 

+  V  (R)-E 

LI  2  dR2  °  J 

|+H(R) 

\|f  =  0  . 


(2.6) 


By  analogy  applying  a  similar  expansion  to  an  electronically  inelastic  process  described  by 
Schroedinger  equation  (2.2)  we  come  to  the  set  of  ordinary  differential  equations  in  the  mixed 
representation: 


dR 


-E 


d  V  . 

Vvv-h  £  jFvvjv'v'W  +  £  U  yyjvv'^)  ^v'v' 

v*v  v  dR  v,s 

(2.7) 


=  0  . 


There  are  two  questions  which  we  need  to  address  here.  First  of  all,  we  construct 
asymptotic  solutions  of  (2.1)  as  h  — »  0  .  It  is  a  relatively  easy  problem  which  can  be  formulated 
in  an  arbitrary  set  of  variables.  It  is  much  more  difficult  to  find  the  connection  formulas  between 
different  asymptotic  solutions.  (As  stressed  by  Schiff11  ^  the  asymptotic  solutions  are  of  little  use 
to  us  unless  we  know  how  to  connect  them  together.)  This  part  of  the  problem  can  be  solved  only 
with  some  additional  assumptions  concerning  the  motions  under  discussions.  Here  we  restrict  our 
discussion  to  systems  with  a  single  unbounded  degree  of  freedom  (a  streamline  motion)  coupled 
with  several  quasi-periodic  degrees  of  freedom  such  that  each  motion  has  a  different  time  scale. 
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There  may  be  also  some  cyclic  motions  which  do  not  create  any  difficulties  on  their  own.  All  the 
bound  motions  are  uncoupled  at  large  values  of  the  streamline  coordinate  R  and  are  assumed  to  be 
adiabadcally  separable  in  the  interaction  region  so  that  one  can  apply  the  usual  WKB  quantization 
rule2. 

Let  us  illustrate  the  made  assumptions  ,  using  scattering  of  a  diatomics  by  an  infinitely 
massive  incorrugated  surface  as  an  example.  The  cyclic  motions  are  represented  by  the  precession 
of  the  diatomic  around  the  normal  drawn  from  the  center  of  mass  of  the  diatomic  to  the  surface. 
Relative  oscillations  of  two  atoms  and  the  nutation  of  the  diatom  give  us  two  quasi-periodic 
degrees  of  freedom .  As  stretching  vibration,  described  by  the  variable  r,  usually  has  a  much  larger 
frequency  than  the  nutation  we  treat  the  former  as  a  quasi-periodic  motion  with  a  frequency 
parametrically  dependent  on  the  nutational  angle  0  and  the  streamline  ccordinate  R.  On  the  other 
hand,  the  nutation  is  expected  to  be  faster  than  streamline  motion  so  that  we  can  describe  it  as 
quasi-periodic  motion  with  the  frequency  dependent  parametrically  on  R.  We  thus  negelect  the 
effect  of  small-amplitude  oscillations  of  the  interatomic  distance  r  on  the  nutation  of  the  diatomic 
and  treat  the  diatomic  as  a  rigid  linear  rotator  in  this  particular  context.  We  shall  see  below  that  one 
needs  these  assumtions  only  to  formulate  the  connection  formulas  for  the  asymptotic  solutions  and 
we  do  not  make  any  approximations  in  the  classical  hamiltonian  itself  to  be  consistent  with  these 
assumptions. 

Note  that  we  included  precession  in  our  analysis  only  to  be  able  to  run  trajectories  in  the 
Cartesian  coordinates.  One  can  directly  start  from  the  Schroedinger  equation  expressed  in  terms  of 
the  curvilinear  coordinates  R,  r,  0,  The  appropriate  kinetic-energy  operator 

a  h2  92 _ h2  92  h2  1  9  .  ft9_ 

2HrdR2  2pr  dr1  2*ir  sin  0  90s,n  90 


is  very  similar  to  that  for  the  J=0  atom-diatom  nonreactive  scattering  problem: 


h2  92 
2^r9R4 


h2  92 

2Hr 


JlV _L  .  J_VJ-  2- 

2  Ur  lirjsinydy 


(2.9) 


with  r  is  the  interatomic  distance  in  the  diatomics,  y  is  the  angle  between  the  diatomics  and  the 
vector  R  drawn  from  the  incident  atom  to  the  center  of  mass  of  the  diatomics  and  R=l  R I.  (Note 

that  the  volume  elements  dR  dr  sin0  d0  and  dR  dr  sinydyare  defined  in  both  cases  in  exacly  the 
same  way.)  Therefore  vibrational  excitation  of  an  diatomics  in  a  J=0  collision  with  an  atom  and 
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scattering  of  diatomic  molecules  by  incorrugated  surfaces  can  be  formally  treated  in  terms  of  the 
nearly-indentical  formalisms. 

To  extend  the  theory  of  semiclassical  transition  probabilites  to  multichannel  equation  (2.1) 
we  start  from  Dirac-Marcus'  representation6,12  of  the  wave  functions  \y^(Q),  ^(Q),  ...,\yn(Q)  as 

Vv(Q)=  Av(Q)  exp[iWv(QVh]  ,  v=l,2 . n,  (2.10) 

where  A  v  (Q)  and  W  y  (Q)  are  real  functions  which  vary  slowly  with  their  arguments. 

Substituting  (2.10)  into  (2.1)  with 

-*(v) 

P  (Q)  =  7  Wy(Q),  (2.11) 


->(v) 

P  (Q) 


2 


=  I  P(V)(Q)  gW  (Q)  P^Q) 

Hit’  p  p 


(2.12) 


2 

and  neglecting  the  terms  of  the  order  of  h  ,  we  come  to  the  following  set  of  differential  equations 
of  the  first  order: 


-E- 


-ihyv 


I  F 
vVv 


-><v') 

p  VV' 

(2.13) 


referred  to  below  as  the  Semiclassical  Multichannel  Equations  (SMEs).  Note  that  div  here 

— > 

implies  use  of  the  covariant  derivatives  of  a  vector  I  so  that 


r«  ig-i/2^gi/2g"i, 

dQ  ^ 


(2.14) 


One  can  easily  verify  that 


1  -Kv)  2  (v) 

IP  (Q)  +Veff(Q)=E 

with 


(2.15) 
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and 


Vevff(Q)HVv(Q)+8Vv(Q) 


(2.16) 


5Vv  S  vV  *  ^  ^ Av  Silt( WV+  Wv)^]/A v  <2‘17> 

and  hence  the  function  W  V(Q)  is  Hamilton's  characteristic  function13  for  effective  potential 

(2.15) . 

The  SMEs  have  the  trivial  solution  for  the  single-channel  problem: 

j\P  (Q)|  +VJ(Q)  =  E  ,  (2.18a) 

r->(D  I 

div\P  p  j J  =  0  (2.18b) 

with 

P  j(Q)  s  A2(Q)  .  (2.19) 

However  what  is  taken  for  free  in  that  case  turns  in  a  rather  complicated  problem  if  one  has 
to  include  interference  between  the  channels.  Even  the  initial  incoming  values  of  effective  potentials 

(2.16)  cannot  be  unambiguosly  determined  because  the  boundary  condition  for  Hamilton's 
principal  function  W  v  (Q)  is  known  only  for  the  initially  populated  channel  labeled  below  by 

index  1.  As  a  result,  quantum-mechanical  correction  (2.17)  to  the  potential  governing  classical 
trajectories  turns  out  to  be  an  ill-defined  function  in  the  infinite- separation  limit. 

One  of  the  ways  to  bypass  this  difficulty  is  to  neglect  quantum-mechanical  correction 

(2.17) ,  compared  with  the  adiabatic  potential  V  V(Q),  bearing  in  mind  that  this  correction  is 

proportional  to  h  and  hence  it  disappears  in  the  limit  h  0 .  In  the  one-dimensional  case  the 
neglect  of  this  correction  leads  us  to  the  set  of  ordinary  differential  equations  that  can  be  solved  by 
means  of  the  finite-difference  method  provided  that  the  coupling  coefficients  are  negligibly  small 
near  the  turning  points  and  hence  one  can  make  use  of  the  standard  WKB  connection  formuale  in 
each  channel.  The  same  assumtion  is  used  by  us  for  the  streamline  motion  in  the  multi-dimensional 
case  however  the  necessity  to  solve  the  set  of  partial  differential  equations  makes  the  problem  much 
more  complicated.  To  overcome  those  difficulties  we  develop  two  approximate  methods  of  its 
solution:  Self-Consistent  Eikonal  Method10  (SCEM)  and  the  multichannel  theory,  Selfconsistent 
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Adiabatically-Corrected  Eikonal  Method  (SCACEM),  both  based  on  using  classical  trajectories  as 

Q 

curvilinear  coordinate  axes  .  In  next  Section  we  outline  this  common  part  of  two  methods,  by 
representing  the  the  SMEs  in  this  particular  curvilinear  coordinate  system  by  means  of  the  standard 
formulas  of  differential  geometry14. 

III.  Traiectorv-following  coordinate  system  (TFCSi 

Let  us  now  express  each  of  Eqs.  (2.13b)  in  terms  of  its  own  set  of  the  curvilinear 
coordinates  q(  v )  =  {t(  v  ),w) with  w  used  for  the  initial  values  of  w  on  trajectories  governed  by  the 

potential  V  (Q)  and  parametrized  by  the  time  t(  v  ^  (instead  of  Q®  =  R)  different  for  each  channel. 

Making  use  of  the  well-known  expression  for  the  covariant  derivative  (see  Eq.(6.87)  in  Ref.  14 
with  g^  for  J  here): 

div  [vv  p  ]%Vv1t7^g(v),v^VVv]  ,  (3.1) 

M  (  v) 


ss' 

where  8(v)  are  the  coefficients  of  momemtum  coupling  in  the  new  set  of  coordinates,  J  v  is  the 
Jacobian  of  the  transformation  from  q^  v)  to  Q  (positive  bv  definition)  and 


awv 

*(v>  • 


(3.2) 


Taking  into  account  that 


with 


.0 

4(V)  ~  ^  * 


(3.3) 


(3.4a) 


for  s  ^  1 


(3.4b) 


we  can  represent  (2.13)  as 
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(3.5) 


ih 


->(v) 

2 

p 

c 


v 


-ihjf  I  F  , 

V  Wv  vv 


->(v) 


•  P 


C  J  J1/2 

V 


where  we  put 


.-,1/2 


Cv(t(v)’w)sJv  (t(v)*w)Vv(t(v)*w)/Al(aw). 


(3.7) 


with  index  1  used  for  the  initially  populated  channel.  Substituting  (2.16),  (2.17)  and  (3.7)  in  (3.6) 
we  find  the  following  equation  for  the  phase  Wy: 


^={E-vf(Q(v,[t(v,.wJ)] 


(3.8) 


and  hence 


wv<*(v) 


l(V) 

,w0)  =  2  /*'{,- v“(Q(vfM.«j} 


(3.9) 


or,  in  a  more  farmiliar  form: 


(v)  H 

Wv<WwO)=  (310) 

Note  that  SMEs  (2.13)  are  independent  of  the  particular  choice  of  the  curvilinear 
coordinates  Q  as  we  assumed  that  the  semiclassical  approximation  can  be  directly  applied  to 
covariant  Laplacian  (2.2).  Use  of  the  Podolsky  transformation15  changing  the  volume  element 
leads  to  a  quantum-mechanical  correction  of  the  order  of  h  which  is  negligible  compared  with  the 
adiabatic  potential  if  the  region,  where  this  correction  is  singular,  is  classically  forbidden.  The 
well-known  Langer  correction16  helps  to  overcome  this  difficulty  for  the  Coulomb  attractive  well. 
In  the  problems  of  our  interest  the  complication  comes  from  the  nutational  motion  because  the 
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appropriate  quantum-mechanical  correction  is  singular  in  linear  configurations,  forcing  the  system 
to  stay  away  from  those. 

To  evaluate  the  Jacobian  Jy  we  make  use  of  the  numerical  algorithm  developed  by  Stodden 

and  Micha17,  extending  it  to  the  equations  of  motion  in  curvilinear  coordinates.  To  be  more 
precise,  we  put 


n^s  _  *(v) 

(v,=  3<I(v) 


ap(v) 

(v)  _  dPn 
- 


(3.11a) 


3qs 


(3.11b) 


(v) 


and  integrate  the  equations: 


.  „ 

Q(V)  = £ 


.  (V)  y 
p  =  I 
ps  ^ 


MM* 


MM 


-  ro  ip(v)  +  r(v)Xo4s8g(v) 

8  [Q(v)]Pp’s  +  pp'^Q(v)9Qp1 


MM' 

,(v)  dg(v) 
p's 


.(V) 


(V)J 

2  MM' 

P's  d  *(V) 


(3.12a) 


v i  Pv  7TO  _ — — 

P'  P  u'W  P'  ^  4’ 

*(v)  *  ^v)**v). 


-IQ 


p's  a 


,  (v)  p  p’ 
*(„)*>(») 
(3.12b) 


together  with  the  equations  of  motion: 


MM 

p(v)  _  T  ag(v)  r(v)  -  ^v- 

P  ~P’^'  P’ 


3Q 


(v) 


9Q 


(3.13) 


(v) 


IV.  Self-Consistant  Eikonal  Method  (SCEM) 


Let  us  now  introduce  the  average  flow  with  the  kinetic  energy  given  by  the  relation 
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1  2  1  -*(v>  2 

2  P  (Q)  =  yIpv(Q)P  (Q)  /p(Q) 


p(Q)SIpv(Q). 


Multiplying  both  sides  of  Eq.(2.15)  by  P  v  ,  summing  over  v  and  taking  into  account  that  the 
— y 

vector  coefficients  F  in  (2.17)  change  their  sign  under  the  interchange  of  v  and  v '  one  can 
easily  verify  that  the  function  P(Q)  satisfies  the  equation 

ip2(Q)+  V(Q)=E  (4.3) 


V(Q)-SV  (Q)p  (Q)/p(Q) 

V  v  v 


and  hence  we  can  define  the  velocity  field  of  the  average  flow  by  means  of  the  relation: 


P  (Q)  =  f7  W  , 


where  W(Q)  is  Hamilton's  characteristic  function  for  the  potential  V(Q) ,  namely. 


W|  +V(Q)  =  E  . 


Let  us  now,  in  following  Micha'  work10,  write  functions  (2.10)  in  the  common  eikonal 


Vv(Q)=Xv<Q>exp[iW(Q)/h]  , 

?v(Q)  =  %v<Q>exrfiW^^’ 


(4.7’) 
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where  we  use  tilda  to  mark  the  diabatic  wavefunctions  by  tilda  to  distinguish  them  from  those  in 
the  adiabatic  represention.  We  neglect  the  difference  between  the  Jacobians  in  Eq.(3.7),  writing 
SMEs  (2.8a),  (3.7)  as 


ih 


dp, 


dt 


(v) 


V  -  V  -  —I 

v  v  2 


►(v) 


Pv  +  I  F  , 

v  V'*v  w 


( -»  -»(v'A 

•  vp+k  ;p. 


V  . 

(4.8) 


where  we  put 


pv(Q(v{‘(vrw; 

D^v^Xv, 

(4.9) 

— Kv)  — Kv) 

— ) 

(4.10) 

o 

III 

(Q)-  P  (Q)  • 

A  similar  equation  in  the  diabatic  representation  takes  the  form 


— >  (v)  2 

K 


2  H  | 

Wv  vv 


V' 


(4.  in 


where  H^,  are  some  diabatic  potential  couplings  and  tilda  is  used  to  distinguish  between  two 
representations. 

The  main  idea  of  the  method  is  to  run  trajectories  on  the  single  potential  surface 
approximating  the  potential  V(Q)  by  means  of  the  relations: 


V(Q)»ZVv(Q) 
v  v 


/P(Q) 


(4.12) 


V(Q) *  I  H  J  (Q)  P*v(Q)/p(Q) 

v  v'  vv  v 

with 


(4.12) 


P(Q)  =  Z  Py(Q 
v1  v 


(4.13) 
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— Kv)  — Kv)  -» 

After  neglecting  the  deviations  K  of  the  momenta  P  from  the  average  P  which  cannot 

be  evaluated  on  the  trajectories  in  question  we  thus  come  to  the  following  equations  for  the 

functions  PV(Q)  and  PV(Q): 


dPv  r  -*  -> 

ihdrJL=Cvv-v]Pv+  v5vFw*  pP. 


i"^=[Hvv-V]|3v+  I  HVV,PV. 


V*V 


(4.14) 


(4.14*) 


Note  that  both  sets  of  equations  conserve  the  right  side  of  (4.13)  along  any  trajectory  governed  by 
the  potential  V(Q) .  The  equations  are  solved  with  the  boundary  conditions: 

Pv(0,w)  =  Pv(0,w)  =  8vlIp0(w)]1/2exp^-ik1R0+  iW°(w)/h)  .  (4.15) 

Compared  with  the  discussion  presented  in  Section  ID,  a  certain  simplification  of  the  formalism 
comes  from  the  fact  that  potential  (4.4)  does  not  contain  the  ratio  of  functions  Ay(Q)  and  \XQ) 

and  hence  one  does  not  have  to  deal  with  the  indeterminate  forms  nccessarvLto  evaluate  the 
effective  potentials  Vy  at  t=0.  It  is  worth  pointing  out  that  the  neglect  of  the  deviations  of  the 

-Kv)  -» 

momenta  P  from  the  average  P  in  Eqs.(4. 11),  (4.11')  is  the  most  serious  assumption  made 
so  far  because  it  changes  the  asymptotic  behavior  of  wave  functions  after  the  couplings  between 
the  channels  completely  turn  off.  In  next  Section  we  show  how  one  can  eliminate  this  deffect  by 
running  a  separate  trajectory  in  the  field  of  each  adiabatic  potential,  instead  of  a  single  trajectory  in 
the  field  of  the  average  potential  V(Q) . 

The  set  of  equations  (4.13')  have  been  recently  derived  by  Stodden  and  Micha18  in  a  more 
tedious.way  by  evaluating  directly  the  time  derivative  of  Jacobian  (3.2).  It  should  be  emphasized 
that  although  similar  sets  of  coupled  differential  equations  have  been  derived  in  the  literature, 
starting  from  the  time  dependent  Schroedinger  equation  in  either  adiabatic19  or  diabatic20 
representations,  the  crucial  difference  comes  from  turning  points  where  the  phases  of  the  functions 
P  V(Q)  and  P  V(Q) :  take  a  discontinuous  change.  The  novelty  of  our  approach  when  extended  to 

the  multi-dimensional  case  comes  from  the  explicit  use  of  the  time-scale  separation  to  carry  those 

functions  through  the  caustics.  We  show  below  that  the  equations  of  motion  used  by  us  in  the 

20 

multi-dimensional  case  also  differ  from  those  derived  by  Meyer  and  Miller  . 
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Potential  (4.12')  has  a  very  interesting  feature,  namely,  its  first  derivative  with  respect  to 
time  is  given  by  the  simple  relation: 

Vs  IP'VH  Jvf'  /p.  (4.16) 

v,v' 

To  prove  it  we  need  just  represent  (4.12*)  and  (4.13')  as 
P 

V  =  P  H  p/p 
and 

ihP=(H-V_l)P 


(4.17) 

(4.18) 


and  then  substitute  (4.18)  in  the  derivative  of  (4.17)  with  respect  to  t,  taking  into  account  that  the 
matrix  H  -  V  1  commutes  with  H  .  The  direct  consequence  of  the  proved  result  is  that  in  the  one¬ 
dimensional  case  the  trajectories  are  governed  by  the  Erenfest  potential: 

.1  h  dH 

P  =  P  - rP/p  .  (4.19) 

-  dQl  " 

(We  can  always  make  g^  *  equal  to  a  constant  1/p^  by  the  appropriate  change  of  variables  so  that 

no  additional  term  appears  in  the  equation  for  the  streamline  motion.) 

The  governing  potential  force  in  the  multi-dimensional  case  has  a  more  complicated  form: 


,P 


P  9H 


8wf 


F  =P  — =n-P/p+  X  n 
-  9Q  ~  P>09Q 


f  P 


LV 


PfM  HP+  P  HP  ^-p^P  HP/p 

(4.20) 


\ 

/P 
/  J 


where  we  put 


Evaluation  of  derivatives  (4.21) 
discussed  in  Appendix  A. 


P  • 


(4.21) 


s 

dw 

along  trajectories  governed  by  Erenfest  potential  (4.12')  is 
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The  most  essential  obstacle  in  solving  Eqs.  (3.6)  in  the  case  of  more-than-one  channels  is 
that  the  functions  £  ,  with  v '  ^  v  in  the  right-hand  side  of  (3.6)  are  supposed  to  be  evaluated 

on  the  trajectory  Q  j.  t  ^  ^ ,  w  ]  .To  overcome  this  difficulty  let  us  first  assume  that  dependence 

— ^  1  2 

of  the  nonadiabatic  couplings  F  ,  on  the  quasi-periodic  coordinates  Q  ,  Q  «...  is  negligible  at 

w' 

large  values  of  the  streamline  coordinate  R=Q®  and  only  the  component  F^v,  along  the 

streamline  motion  is  asymtotically  important.  As  a  result  of  this  assumtion  SMEs.(2.13)  become 
asymptotically  separable.  This  implies  that 


>1  .. 

JR=DO 


hkv(w>ar 


(where  k,  is  just  a  constant)  and  hence 


arg[¥v(°-w)]  =  ~  ikv(w)Ro+  +  mod(rc)  . 


P^)(w)  =  P((l!)(w)  H  >0  , 


:v(w)  =^2(lR  [V^RofW)  -  vv(Rorw)]  +  h  kl  . 


We  represent  now  Eqs.(3.6)  as 


ih't(v-)[R’w)('5Tr)w=[Vv( 


x  1 1 — >(v— )  2 

Q(v  _)[R,WV  “  E  -  yl  p  IR,W]  JCytR.w] 


J^2[R,w]  -4  /  x  ->(v)  i5WuTQfv  ,1/h 

■ihvw7^~ ;fwK-)(R-w))*p  [Q(VJ? 

J^Q(V->]  ,5.5, 
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where  minus  marks  the  incoming  wave,  the  components  P^  [Q^v  of  the  momentum 
—K  v) 

P  [Q(V  _)]*  Jacobian  J V'_[Q(V _)]  an<*  the  phase  shift  8Wv,[Q(v_)]  are  approximated 
respectively  as 


'  Ou.  (v) 

8  /2  p0  [Q(v.)]=- P(v, 


(5.6a) 


v  [Q(v-)]-p; 


( vl 

-Cevj 


(5.6b) 


±[Je  -  v(R,Q>...,QMQ;y+1,...)]-7E  -  v(R.Qt....,QM.QPv),QP+I....) 


„  .  1  y  <¥p<v) 

P<V>  2nto*  I*  ’ 


U.'  1  (  u'  p'  ^ 

Q  [  R,  w]  ■  y|Q(  y  _}( R,  w]  +  Q^.  }[  R,  w]  J, 


P[tV’  [Q(v-)1 

VtV-)] 


v  ;  U  (v) 

"'vWvj]-^,  l  d<5  P4  [Q(V_)3  • 

Q(V) 


(5.10) 
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(Sign  (+)  or  (-)  in  (5.6b)  is  chosen  depending  on  the  direction  of  the  appropriate  one-dimensional 
motion.)  We  integrate  (5.5)  by  running  a  set  of  trajectories  on  the  adiabatic  potential  surfaces  with 
the  initial  conditions  given  by  Eqs.(5.3),  (5.4)  for  the  momenta  and 

Cv((Xw) =  5vj[p°(w)]1^2exp(-  UcjRq  +  iW°(w)Ai)  n 

for  the  wavefunctions  themselves. 

The  most  important  consequence  followed  from  the  presented  mathematical  arguments  is 
that  the  semiclassical  multichannel  solution  sought  for  does  have  the  eikonal  form  —  remember  that 
existence  theorems  have  been  proved  only  for  single-channel  processes12  .  The  equations  for  the 
outcoming  wave  are  similar  to  (5.5)  except  that  they  are  solved  with  a  diifferent  initial  conditions 
formulated  using  the  WKB  connection  formula  at  the  classical  turning  point  for  the  streamline 
motion  on  the  appropriate  adiabatic  potential  surface. 

VI.  Evaluation  of  the  S-matrix  elements  from  semiclassical  wavefunctions 

To  avoid  the  singularities  at  the  classical  turning  points,  when  evaluating  the  S-matrix 
elements: 


S 


(v) 

v,v 


f  I 

1)=±1 


JdwVvv(Rf w,v)uv(w)  . 


(6.1) 


where  u^/w)  is  an  exact  quantum-mechanical  wave  function  and  v  is  a  set  of  the  appropriate 
quantum  numbers,  we  first  represent  (6. 1)  as 


s<v> 


1 

AR 


Rf+AR 


JdR  S 


(v) 

v;v' 


(6.2) 


The  relation  turns  into  the  identity  if  we  deal  with  the  exact  matrix  element  which  is  independent  of 
R.  We  thus  find 


AR 


(v)  (v) 
v  Sv;v'" 
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r-r-r -  1  f(w0-'l))  jj2 

„?±1  liXJk  <w0-',)  l  d‘  Jv  ('•wO'’,)?’w<,'wO’1))Uv'(*(1);'vo)  ’ 

0  (6.3) 


where 


C,w(t’wav)sCw(t’w0>'))e 


-^L  =ps(ws;ws  t))/m 

dr  u 


-ik(v)(w0,i))R(t,w0,u) 


(6.4) 


(6.5) 


and  J  v  is  the  absolute  value  of  determinant  of  the  auxiliary  matrix  with  the  elements 


u) 


propagated  instead  of  the  functions 


(6.6) 


Us 

to  avoid  singularities  in  the  initial  conditions  for 


functions  (3.11b).  (The  function  P^(w)  in  (5.11)  is  chosen  in  such  a  way  that  the  matrix  element 

(1) 

S  v.  v  is  equal  to  1  for  elastic  collisions.)  The  crucial  advantage  of  the  integral  (6.3)  compared  with 
(6.1)  is  that  the  Jacobian  J  v  appears  now  in  numerator  whereas  function  (6.4)  is  nonsingular. 


CQnclusigflS 


The  following  formal  issues  were  examined: 

(i)  development  of  prescriptions  for  continuation  of  global  semiclassical  wavefunctions 
through  classical  turning  points, 

(ii)  development  of  a  practical  numerical  algorithm  to  extract  final  populations  of  vibrational 
levels  from  the  semiclassical  outgoing  solution, 

(iii)  exploration  of  means  to  go  beyond  the  assumptions  underlying  SCEM  to  take  into 
account  specific  topology  of  the  adiabatic  potential  surfaces,  especially  responsible  for  tunneling 
and  related  resonances  which  gave  rise  to  the  new  multichannel  semiclassical  theory  referred  to  as 
SCACEM. 
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Detailed  numerical  studies  such  as  the  ones21  in  progress  are  expected  to  reveal  the 
potential  and  power  of  this  framework. 

Appendix  A 

The  derivatives  of  w  ^  with  respect  to  Q  are  elements  of  the  matrix  Q'^with  Q  formed 

sss  « 

by  the  functions 


(Al) 


We  thus  come  to  a  nonlinear  set  of  differential  equations  —  we  need  to  know  the  derivatives  of  Q 
with  respect  to  q  to  write  the  equations  of  motion  and  we  need  the  equations  of  motion  to  evaluate 
those  derivatives.  The  situation  is  complitely  different  from  that  for  the  single-channel  problem 
because  in  the  latter  case  the  potential  is  the  known  function  of  coordinates.  To  decouple  the 
equations  we  neglect  the  second  derivatives  of  the  functions  P  ,  P  and  Q  with  respect  to  q.  One 

finds 


.(v)  \i\i  n 

p..  =xq,:j 


.  t-  a  h 


P  ~£'<(v)K  n  n" 

P  dQ  dQ 


(  ^ 

+  2Rq  P 


aq 


PK,P 


P  dH 

VP-P.U  P  “^h-P/P 


J 


'  ,IX  K  IX 

^  -  a<f 


+  2  I 


dw 


p' 


P^dQ1 
P 


P 


Re<P,p  5P,^)_p^'Re(P  gP,p>0 


/P+P^P.vP  jP/p 


_3 


+  X 


dw' 


P’>0  d  Q4 


T  ••••" 

r  1-  dH  fan  1“  d  H  ' 

k  £ 

a 

— TT - 1 

P  p  irP  +  P  ~=\XT  P  Li  P,pP  LX*  P/p 

r’  aq  “  ”  aQ  “  aq  J 

/p 

(A2) 


with 


pW.’ 

<v) 


dq 


(A3) 


(v) 
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Functions  (Al)  are  then  found  by  integrating  the  equations: 


.  W*’  wt’  • 

6  =  p 

^(V)  r(v) 


(A4) 


To  take  into  account  the  divergence  of  streamlines  we  use  Eqs.(4.16)  only  to  find  the  phase 
of  the  function  %y(Q)  whereas  its  absolute  value  |PV(Q)|  is  evaluated  by  solving  the  equations: 


£(,Pv)=j  ^vta(xv,xi) 


<A5) 


with 

Ev^flV'Pv1 


and  J  used  for  determinant  of  the  matrix  Q: 

s 

"  a(q)  ■ 


(A6) 


(A7) 


A  similar  scheme  can  be  developed  in  the  adiabatic  representation. 
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ACCURATE  COMPUTATIONS  OF  ENERGY  TRANSFER  IN  ATOM-DIATOM 
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(PRELIMINARY  DRAFT  for  Journal  of  Chemical  Physics) 


by 

P.  K.  Swaminathan,  G.  A.  Natanson,  B.  C.  Garrett,  and  M.  J.  Redmon 
Chemical  Dynamics  Corporation 
Upper  Marlboro,  MD  20772 


I.  INTRODUCTION 


Molecular  collisional  energy  transfer  presents  a  computationally  challenging 
problem  to  the  dynamicist  even  after  the  global  potential  energy  surface  is  known.  A  rich 
variety  of  quantum  mechanical,  semiclassical  and  classical  models  may  be  constructed  to 
solve  the  problem  depending  on  the  behavior  of  the  potential  energy  surface  which 
controls  the  collision  dynamics.  The  information  from  modem  laser-based  experiments, 
involving  precise  microscopically  detailed  observations  of  collisions,  is  in  the  form  of  state- 
resolved  cross  sections  and  demand  an  explanation  from  theory  at  that  level  of  detail. 
Classical  mechanics  in  itself  does  not  offer  a  long-term  solution  for  obtaining  such  state- 
resolved  information.  Fully  quantum  mechanical  treatments  become  unwieldy  for  realistic 
systems  despite  recent  triumphs  that  exploit  efficient  algorithms  on  supercomputers.  Thus 
there  is  a  present-day  need  to  develop  and  exploit  semiclassical  dynamical  theories  that  hold 
promise.  In  order  to  obtain  general  results,  it  has  to  be  based  on  a  systematic  process 
involving  detailed  benchmarks  against  accurate  quantum  mechanical  results  for  smaller 
systems  and  the  development  of  controlled  computational  schemes  which  arc  careful  paths 
to  accomplish  the  devious  transcription  from  classical  trajectories  to  semiclassical 
wavefunctions. 

The  present  paper  reports  a  study  of  vibrational  excitation  during  atom-diatom 
collisions  in  gas  phase  and  is  one  of  many  successful  numerical  experiments  in  which  we 
have  employed  the  eikonal  ansatz  for  the  semiclassical  wavefunction.  The  relevant  general 
formal  discussions  have  been  presented  elsewhere  in  greater  detail.  The  present  article  will 
focus  on  the  formulation  for  the  vibrational  excitation  problem  and  its  numerical 
investigation. 

In  this  paper  we  consider  a  two  degree  of  freedom  model  of  atom-diatom  collisions 
(O+HF  and  Na+H2  collisions  are  studied)  as  defined  in  the  rotational  infinite  order  sudden 

(IOS)  approximation.  The  numerical  tests  presented  here  are  benchmarked  against 
previously  existing  quantum  IOS  (QIOS)  results.  Two  distinct  dynamical  models  are 
employed  in  this  paper  to  describe  state-to-state  vibrational  dynamics:  (1)  is  based  on 
combining  a  classical  treatment  of  translational  motion  with  a  quantal  states  expansion  in 
target  vibrational  states  (called  RSVE:  for  rotational  sudden  with  vibronic  expansion)  and 
(2)  is  based  on  employing  semiclassical  eikonal  wavefunctions  for  all  degrees  of  freedom 
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obtained  from  their  classical  trajectories  (called  SIOS:  for  semiclassical  infinite  order 
sudden). 

n.  THEORY 

l.IOS  Approximation  in  the  diabatic  representation 
The  rotational  IOS  approximation1  for  the  atom-diatom  scattering  problem  is  based 
on  solving  the  two-dimensional  Schroedinger  equation  with  the  Hamiltonian: 


A 

h 

IOS 


hV 


2  \i 


dRJ 


fa2  a2  h2i(i  +  i) 

2jiBC  9r2  2pR2 


fa  2  j(  j  +  1)> 
2^BCr2  > 


I +  H(R,r,y) 


(1.1) 


where  H(R,r,y)  is  the  n  x  n  matrix  of  the  diabatic  potentials  and  couplings  between 

different  electronic  states,  R  is  the  distance  between  the  center  of  mass  of  the  diatomic  BC 
and  the  atom  A,  r  is  the  bond  length  in  the  diatomic  BC  and  7  is  the  angle  between  the 
diatomic  and  the  radius  drown  from  its  center  of  mass  to  the  atom  A. 


2.  Quantum-mechanical  treatment  of  vibrations. 

In  following  Redmon  et.  al.  1  we  reduce  the  problem  to  the  infinite  set  of  the 
ordinary  differential  equations  of  the  second  order  by  using  the  matrix  representation  for 
the  hamiltonian  in  the  space  of  the  vibrational  coordinate  r,  namely,  we  use  the 
eigenfunctions  of  the  Schroedinger  equation: 


2  2 

h  d 


^rjT  +  VBC(r)  + 


h2Xj  +  1)  j 

SUbc'2  £v. 


ui  =  0 


as  the  basis  set  and  represent  (2.1)  as 


+  X  h 
n’.v' 


jl 

nv.n’v 


=  0  , 


(2.1) 


(2.2) 


2 


where  we  put 


jl  1(1+ 1) 

h  ,,(R,Y)=S  ,8  ,  --0— 

nv,nV  '  nn  w  R2 


+  7t{-  ADTl+^dr  uJv(r)uJv-(r)[Hnn^r*Y)  -SnnlVBC(r)]}  ’ 

ti 

(2.3) 


with 


hkjv  *  [2(i(E  -  eJv)] 


1/2 


(2.4) 


and 


AD„  = 


(2.5) 


used  to  evaluate  the  probabilities  of  transitions  from  the  T)th  electronic  state. 

We  solve  the  set  of  coupled  differential  equations  (2.2)  by  representing 
wavefunctions  ^v  in  the  common  eikonal  form3: 


v^v(r;y)  = 


1 


-Pnv<R;Y>cxP[W  (R;YVh] 


p  (R;y) 


(2.6) 


where  (R;y)  is  Hamilton's  characteristic  function4  found  from  the  Hamilton-Jacobi 
4 

equation  : 


3 


(2.7) 


1  I  jl  p  _jl  12 

2TTlP  <*y)|  +  v  (R;y)  = 


with 


Jl_dW 
dR 


jl 


(2.8) 


_jl 

and  V  (R;y)  denotes  the  Ehrenfest  potential: 


Vjl<R:r)  =  „v^v'-.nVp!vPnv‘/S 


P' 


nv 


(2.9) 


jl 

The  functions  Pny  satisfy  the  set  of  coupled  equations: 


jl  fl  fl  jl  a 

ih^V%vV"v.nV'PnV-V  P"V 


(2.10) 


and  are  solved  together  with  the  equation  of  motion  for  the  potential  (2.9). 


3.Semiclassical  treatment  of  vibrations 
A.Equations  of  Motions 

We  represent  the  wavefunction  V(Rn  Ylj)  in  the  eikonal  form: 
\jf(R,r,r,lj)  =  A(R,r;y,lj)  exp[iW(R,r;y;lj)/h] 

where  W(R,r;y;lj)  is  Hamilton's  characteristic  function: 
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(3.2) 


for  the  potential 


U(R,„.j,.V(R,r:r)  +  iL!^  +  i^ 

2t*R  2HBCr 


The  momenta  conjugate  to  the  coordinates  R  and  r  are  given  by  the  relation: 


P0(R,r;Y;S)  =  ^-,P1(R,r;r.S)=  5^. 


To  cany  the  wavefunction  V(R»  r,  y, lj  )  through  turning  points  we  represent  it  as 
-1/2 

V(t,r0;y,lj )  s  J  (t,r0;r,lj )  A(ftr  ;j)  C(t,r0;y;lj ) 

where  t  is  time,  Tq  is  the  initial  value  of  r  on  the  trajectory  run  on  the  potential  surface  U 
and  J  (t,r^;  y,lj )  is  the  Jacobian  of  the  transformation: 

t,r0  ->  R(t,rQ;y;ljX  r(t,rQ; y;lj)  ,  (3.6) 


and  then  assume5  that  the  function  ^(t.r^iyjlj )  is  continuous  everywhere  by  analogy  with 
the  one-dimensional  case.  We,  in  following  Ref.5,  solve  the  equation: 
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(3.7) 


as  if  the  function  (^t.r^yjlj )  itself  were  continious,  finally  shifting  its  phase  by  ji/4  times  the  total 
number  of  the  turning  points  passed. 


(3.8) 

(3.9) 

(3.10) 


with 


Qj(0,ro,±  1)50,  Ql(0,ro»  ±  1)  H  1 


(3.11) 

(3.12) 

(3.13) 
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P®(tu(r0,±  1)  =  0,  P^O.r^l) 

•*-n 

% 

=L 

1 

1 

P1(0,r(),±  1)  dr| 

r=rn  , 

0  (3.14) 

r0  1 
W0(r0,±  l)  =  ±^/dr'pr(r*)-f 

* 

(3.15) 

f(0,r0,±  l)s-p^./p  , 

(3.16) 

WV*  1)=4[kyjlPr(r0)lf1/2 

iW0(r0,±l)-ik  R0 
e 

(3.17) 

i  tn  iW„(rA,±  l)-ik  .Rn 

C(0,r0,±  l)  =  ^IJ(ar(y±  l)l1/2e  VJ  °. 

(3.18) 

c.Connection  formulae 

If  the  interaction  between  vibrational  and  collisional  motions  is  negligible  the  derivative  of  r 
with  respect  of  its  initial  value  r^  can  be  evaluated  by  differentiating  the  equation: 

r(t,r0) 

t=,lBC  J  dr7pr(r'>  (3.19) 

0 


with  respect  to  Tq  at  the  fixed  moment  L  We  find 


(3.20) 


and  hence 
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Ol‘^jt=pr(r>/Prfr0) 


(3.21) 


Therefore  if  the  off-diagonal  element  can  be  also  neglected  the  Jacobian  J  factors  as 

HJ(t,r0,±  l)»P°(t,ra±  l)P1(t,r(y±  1)/P1(t()fr(y±  1)  ,  (3.22) 

and  its  zeros  approximately  coincide  with  zeros  of  the  momenta  P®  and  P*.  We  thus  assume  that 
jumps  of  the  Jacobian  on  the  caustics  are  governed  for  adiabatically  separable  motions  by  the 
standard  WKB  rules,  namely,  that  the  function  l£(Q)l  is  continuous  there,  whereas  the  phase  of  the 
function  C(Q)  decreases  by  n/2: 

W(Q)->W(Q)-a/2  (3.23) 


if  the  Jacobian  changes  its  sign,  and  by  n,  if  both  momenta  change  their  sign  during  the  same  time 
step. 


The  exact  quantum  mechanical  wave  function  has  the  asymptotics: 


"vj<r)  -ikvjRf 

n/Rrr>  =  Ti73 + 


vj;v5' 


"vi(r>  *.TRf 


(3.24) 


where  the  S-matrix  elements 


for  vj  *  v'j'  are  given  by  the  relation: 
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By  analogy  with  (3.25)  we  evaluate  the  semiclassical  S-matrix  elements  by  means  of  the  relation: 


i  _ _  -ik.,Rf  — 

S  .  ,,=  ~-  / k  e  1  £  /dry  .(R„r;\))u  ,.(r)  .  (3.26) 

vjjv'j  2-y /  v’J  v=±l  q  vjv  f  vj 

To  avoid  of  the  singularities  in  the  turning  points  we  first  represent  this  matrix  element  as  the  two- 
dimensional  integral  using  the  relation: 


Rf+AR 


S  .  JdR  S  . 

vj;v'j’  AR  1  yj;vj 


(3.27) 


(The  relation  turns  into  the  identity  if  we  deal  with  the  exact  matrix  element  which  is  independent  of 
R.)  Substituting  (3.26)  into  (3.27)  we  find 


1/2, 


AR/k  A  S  . 

V  vj^BC  yj )  vj;v  j 


(3.28) 


J'±lrJ-dr°Lpr(r(Tu) 


r  k(r  u)  *f(ro»')) 

°  ^  J1/2Ct,r0>t>)C,yj(t,r0,v)uv3.(r(t);r0) 


T/2  1/2 

=  -^4  J  dxk^Uu]  /  dtT  [t,x,\>] C"yj[t,'i;v]uv^,(r(t);r0(T)) 


v=±l  0 


0 


where 


and  the  notations  tx.u]  and  [t,t,\)]  stand  for  (rg(t),\))  and  (t,  r^Cx),!)).  We  directly 
0  1 

propagate  Qj  and  Qj  along  the  trajectory  to  avoid  possible  complications  coming  from  the 
singularities  in  ( 3.14)  which  appear  if  the  trajectory  starts  from  the  turning  point 

III.  POTENTIAL  ENERGY  SURFACES 

Two  candidates  of  collisional  excitation  of  vibrational  levels,  (1)  in 
O(3P)+HF(v=0,j=0)  collisions  and  (2)  in  Na(3s2S)+H2(v=0,j=0)  collisions,  both  of 

which  have  previously  existing  potential  energy  surface  information,  were  studied.  The 
potential  energy  surface  for  the  OHF  system  is  based  on  an  earlier  UHF  calculation2  and  a 
sum  over  pairs  plus  three-body  correcting  analytical  fit  (surface  1  of  ref.  2)  whereas  the 
NaH2  collisions  were  studied  on  a  DIM  potential  energy  surface.6  A  LEPS  fit  to  MCSCF 

results  for  OHF  were  also  generated  in  ref.2  to  enable  the  dynamical  study  of  the  H 
abstraction  reaction  but  this  needs  to  be  employed  only  at  higher  energies  than  presently 
reported.  We  refrain  from  including  the  analytical  and  numerical  details  regarding  the 
potential  surfaces  for  brevity  since  they  are  published  elsewhere. 

A  central  feature  of  the  OHF  potential  surface  that  is  relevant  for  vibrational 
excitation  dynamics  is  the  expected  concentration  of  most  of  the  vibrational  excitation  (at 
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lower  energies)  to  a  narrow  range  of  IOS-angles,  close  to  the  approach  of  O  atom  from  the 
H-atom  end,  of  the  HF  molecule;  in  this  direction,  the  potential  surface  allows  a  reactive 
encounter  to  occur  at  higher  energies.  This  will  be  borne  out  by  the  dynamics  results 
below.  In  the  case  of  NaH2  collisions,  only  the  electronically  adiabatic  ground  surface 

(obtained  from  diagonalizing  the  8  diabatic  DIM  surfaces  of  ref.  6)  is  employed  in  the 
electronically  elastic  results  to  be  presented  below  for  comparing  QIOS  and  RSVE  results. 

IV.  CALCULATIONAL  DETAILS 

The  RSVE  dynamical  calculations  using  eqs.  (2. 10)  are  straightforward  once  the 
potential  matrix  in  the  vibronic  basis  is  setup  to  be  computed.  This  is  done  by 
diagonalizing  the  asymptotic  vibrational  Schroedinger  equation  in  a  harmonic  basis  set  or 
employing  one  of  many  standard  numerical  algorithms1  to  obtain  the  target  states.  The 
atom-diatom  potential  matrix  is  readily  computed  along  the  trajectory  since  it  is  expanded 
in  these  target  states  which  have  amplitudes  that  vary  with  time  according  to  eq.  (2.10) 
thereby  leading  to  a  time-dependent  potential  obeying  (2.9).  The  amplitude  initial 
conditions  are  chosen  to  correspond  to  unit  probability  in  the  diatomic  v=0  state.  The 
separation  R  is  chosen  large  and  the  asymptotic  negative  relative  momentum  is  defined  by 
the  channel  potential  in  the  v=0  channel  and  collision  energy.  The  final  amplitudes  in  other 
state  channels  resulting  from  solving  the  dynamical  equations  yield  the  inelastic  transition 
probabilities.  The  computation  of  cross  sections  involves  integration  over  many  trajectories 
that  represent  varying  IOS-angles  and  impact  parameters  for  each  given  initial  rovibrational 
state.  All  the  preliminary  computations  reported  here  are  for  initial  v=j=0  for  the  diatomic 
and  the  study  examines  dependence  on  collision  energy  and  total  angular  momentum 
(semiclassically  defined  by  the  impact  parameter). 

The  SIOS  calculations  start  with  the  generation  of  a  bundle  of  classical  trajectories 
determined  by  the  collision  energy  and  total  angular  momentum.  Only  20  trajectories  are 
employed  in  the  present  results.  The  main  task  involves  building  the  semiclassical  eikonal 
wavefunction  along  these  trajectories  giving  special  attention  to  the  semiclassical  phase 
changes  at  caustics  (such  as  classical  turning  points).  The  latter  are  determined  by  the 
semiclassical  Bom  Oppenheimer  approximation  (see  previous  section  and  Appendix  2) 
applied  to  the  present  problem.  The  post-collision  eikonal  wavefunction  is  then  projected 
onto  selected  final  states  to  directly  determine  the  S-matrix  elements  for  inelastic  transitions. 
The  projection  process  employs  interpolation  of  the  asymptotic  eikonal  wavefunction  to 
result  in  a  finely  grided  numerical  quadrature  for  S-matrix  amplitudes  since  the  inelastic 
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trajectories  end  in  a  nonuniform  grid  even  when  the  initial  conditions  start  from  a  uniform 
grid.  Cross  section  computations  involve  integration  of  these  amplitudes  over  IOS-angles 
and  impact  parameters  and  introduce  further  quantal  interference  effects. 

V.  RESULTS  AND  DISCUSSION 

Figure  1  shows  that  excellent  agreement  with  an  exact  limit  QIOS  benchmark  is 
obtainable  from  the  RSVE  computations.  This  success  of  RS  VE  in  gas  phase  is  consistent 
with  the  similar  success  that  the  coupled  states  ansatz  has  had  in  treating  rotationally- 
electronically  inelastic  gas/surface  encounters  involving  up  to  282  channels.7  The  most 
exciting  aspect  of  figure  1  is  however  the  promise  shown  there  by  the  SIOS  method,  which 
is  seen  to  capture  the  essential  physics  of  this  collision;  the  results  shown  are  from  SIOS 
calculations  that  employ  semiclassical  wavefunctions  based  on  a  meager  20  trajectories  and 
are  not  yet  numerically  the  best  they  can  be. 

Figures  2  (a  and  b)  show  dependence  on  total  angular  momentum  for  the  two 
examples  of  vibrational  excitation  studied  here.  Fig.  2a  contains  all  the  data  available  at  the 
time  of  this  report.  (Detailed  calculations  are  in  progress  to  make  a  full  assessment  of  the 
methodologies  in  different  L  and  E^  regimes.)  The  NaH2  computations  were  done  to  see  if 

the  quality  of  agreement  is  maintained  for  a  completely  different  molecular  system  and  Fig. 
2b  is  encouraging  (note  that  H2  being  a  homonuclear  diatomic  molecule  is  fully  studied  by 

sampling  half  the  IOS-angle  space  compared  to  HF). 

VI.  CONCLUDING  REMARKS 

We  have  produced  a  systematic  numerical  benchmark  for  collisional  excitation  of 
vibrational  levels  employing  semiclassical  methods.  Two  semiclassical  computational 
routes,  (1)  using  a  coupled  quantal  states  ansatz  to  describe  vibrational  states  quantum 
mechanically  and  translational  motion  classically,  within  the  SCEM  framework,  and  (2) 
using  semiclassical  wavefunctions  generated  from  classical  trajectory  results  using  a  newly 
developed  (semiclassical  Bom-Oppenheimer  approximation)  procedure  were  tested  against 
exact  quantum  limit  results.  Both  the  semiclassical  methods  were  successful  in  obtaining 
reliable  numbers  for  state-to-state  vibrational  excitation  probabilities.  The  use  of  method 
(2),  termed  SIOS  herein,  is  indicated  as  an  exciting  prospect:  such  a  technique  is  to  be 
welcomed  as  a  computational  route  allowing  the  transformation  of  information  from 
classical  trajectories  to  state-to-state  transition  amplitudes  at  the  semiclassical  level.  Tiie 
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novel  feature  of  the  SIOS  scheme  employed  here  compared  to  rigorous  S-matrix  theory  of 
Miller®  (which  was  difficult  to  apply  in  practical  problems  and  soon  went  out  of  use)  is 
that  there  are  no  nonlinear  root  trajectory  searches  to  be  performed  here;  instead,  amplitudes 
are  obtained  by  projecting  the  eikonai  scattering  wavefunction  on  to  specific  final  states. 
The  promise  of  such  a  method  was  demonstrated  earlier  in  the  Franck-Condon  region  for 
photodissociation,9  but  the  present  validation  for  vibrational  excitation  processes 
significantly  extends  the  validity  regime,  by  demonstrating  that  useful  quality  persists  in 
outgoing  eikonai  wavefunctions  in  the  asymptotic  region  provided  valid  connection 
formulae  are  employed. 

ACKNOWLEDGEMENT 

Support  of  this  work  from  AFOSR/SBIR  contract  #F49620-88-C-0085,  is 
greatfully  acknowledged. 


13 


REFERENCES 


1  See  eg.  V.  Khare,  D.  J.  Kouri,  J.  Jellinek,  and  M.  Baer,  in:  Potential  Energy 
Surfaces  and  Dynamics  Calculations,  D.  C.  Truhlar,  rd.,  (Plenum,  New  York,  1981) 
p.  475. 

2  M.  J.  Redmon,  L.  T.  Redmon,  and  B.  C.  Garrett,  in  AFRPL  Final  Report  AFRPL 
TR-84-030,  "Collisional  Excitation  Cross  Sections" ,  Sec.  3.  (1984). 

3  D.  A.  Micha,  J.  Chem.  Phys.  78,  7138  (1983). 

4  H.  Goldstein,  "Classical  Mechanics",  Addison- Wesley,  London  (1950). 

5  G.  Natanson  and  P.  K.  Swaminathan,  "Semiclassical  Methodology  for  Computing 
Multichannel  Eikonal  Wavefunctions  in  Molecular  Collisions:  A  Reformulation  and 
Extension",  to  be  published.  A  draft  of  this  is  Appendix  2  of  this  Phase  I  report 

6  D.  G.  Truhlar,  J.  W.  Duff,  N.  C.  Blais,  J.  C.  Tully,  and  B.  C.  Garrett,  J.  Chem. 
Phys.  77,  764  (1982);  N.  C.  Blais,  D.  G.  Truhlar  and  B.  C.  Garrett,  J.  Chem. 
Phys.  78,  2956  (1983). 

7  B.  M.  Rice,  B.  C.  Garrett,  P.  K.  Swaminathan,  and  M.  H.  Alexander,  J.  Chem. 
Phys.  90,  575  (1989). 

8  See  W.  H.  Miller,  Adv.  Chem.  Phys.  25, 69  (1974)  for  a  review. 

9  P.  K.  Swaminathan,  C.  L.  Stodden,  and  D.  A.  Micha,  J.  Chem.  Phys.  In  Press 
(1989). 


14 


Figure  Captions 

1 .  Comparison  of  angle-dependent  probabilities  for  v=0  ->  v=l  vibrational  excitation 
in  O(3P)+HF(v=0  j=0)  collisions  at  E^S.OeV,  for  total  angular  momentum,  L=0. 

Points  arc  marked  by  bullets  =  SIOS,  circles  =  QIOS,  and  stars  =  RSVE  results. 

2.  (a)  Comparison  of  RSVE  and  QIOS  results  for  variation  of  IOS-angle  dependent 
v=0  ->  v=l  probabilities  with  total  angular  momentum  L  at  Etr=3.0eV  for 

O(3P)+HF(v=0,j=0)  collisions. 

(b)  Comparison  of  RSVE  and  QIOS  methods  for  variation  of  IOS-angle  dependent 
v=0  ->  v=l  probabilities  with  total  angular  momentum  L  at  Etot  =  0.075  au  for 

Na(3s2S)+H2(v=0  j=0)  collisions. 
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